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ABSTRACT 

A  new  remote  sensing  approach  based  on  polarimetric  wavelet  fractal  detection  principles  is  introduced  and  the  Mueller 
matrix  formalism  is  defined,  aimed  at  enhancing  the  detection,  identification,  characterization,  and  discrimination  of 
unresolved  space  objects  at  different  aspect  angles.  The  design  principles  of  a  multifunctional  liquid  crystal  monostatic 
polarimetric  ladar  are  introduced  and  related  to  operating  conditions  and  system  performance  metrics.  Backscattered 
polarimetric  signal  contributions  from  different  space  materials  were  detected  using  a  laboratory  ladar  testbed,  and  then 
analyzed  using  techniques  based  on  wavelets  and  fractals.  The  depolarization,  diattenuation,  and  retardance  of  the 
materials  were  estimated  using  Mueller  matrix  decomposition  for  different  aspect  angles.  The  outcome  of  this  study 
indicates  that  polarimetric  fractal  wavelet  principles  may  enhance  the  capabilities  of  the  ladar  to  provide  characterization 
and  discrimination  of  unresolved  space  objects. 

Keywords:  Space  Surveillance,  Ladar  Remote  Sensing,  Physical  Algorithm,  Space  Materials  Characterization,  Mueller 
Matrix  Decomposition,  Polarimetric  Wavelet  Fractal  Analysis, 


1.  INTRODUCTION 

Polarimetric  sensing  and  imaging  offer  unique  advantages  for  a  wide  range  of  detection  and  classification  problems  due 
to  the  intrinsic  potential  for  high  contrast  in  different  polarization  components  of  the  backscattered  light.  Indeed, 
polarized  imaging  can  yield  high-specificity  images  under  high-dynamic  range  and  extreme  condition  scenarios,  in 
scattering  media,  or  cluttered  environments,  offering  at  the  same  instance  information  related  to  the  object  material 
composition  and  its  surface  characteristics  1"11,18"30.  Active  polarimetric  sensors  using  laser  sources  offer  predictable 
illumination  levels,  polarimetric  filtering  to  minimize  specular  reflections,  and  the  ability  to  operate  independent  of 
sunlight  illumination  and  temperature.  In  addition,  an  active  illumination  source  is  aligned  with  the  detector,  allowing 
more  direct  light  to  reach  the  objects,  which  in  turn  reduces  shadows  and  increases  visibility,  reducing  at  the  same  time 
the  false-alarm  rate.  Although  active  sensors  offer  superior  signal-to-noise  ratio  (SNR)  and  convenience  of  use,  they  are 
larger,  more  costly  and  complex  to  build,  and  consume  more  power.  While  there  exist  a  number  of  effective  space 
surveillance  detection  and  imaging  spectral-polarimetric  technologies  and  methodologies,  further  progress  is  both 
necessary  and  possible  toward  more  efficient  and  reliable  detection,  identification,  characterization,  and  discrimination 
of  small,  low-contrast,  unresolved  space  objects.  Accurate  sensing  of  small  and  distant  satellites  as  well  as  the 
discrimination  of  these  satellites  from  natural  near-earth  objects  and  man-made  debris  is  still  a  challenging  and  important 
space  remote-sensing  task,  along  with  the  prediction  of  interactions  between  and  threats  to  these  space  objects.  New 
remote  sensing  detection  principles  based  on  polarimetric  wavelet  fractal  detection  principles  aimed  at  enhancing  the 
detection,  identification,  characterization,  and  discrimination  of  unresolved  space  objects  at  different  aspect  angles  were 
previously  introduced  9"10,28.  The  combination  of  polarimetry  and  wavelet- fractal  analysis  yields  enhanced  knowledge  of 
the  spatial-temporal-frequency  characteristics  of  backscattered  ladar  signals  from  space  objects  28. 

Wavelets  are  a  family  of  transforms  whose  basis  functions  are  of  short  duration  and  finite  energy  10.  There  is  a  large 
variety  of  wavelet  functions  available.  Various  wavelet  families  exhibit  different  degree  of  spatial  compactness  and 
smoothness.  In  contrast  to  the  Fourier  transform,  which  effectively  assumes  that  the  signal  is  stationary  at  time  scales  of 
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interest,  the  wavelet  transform  determines  a  signal’s  frequency  content  as  a  function  of  time.  The  result  is  a  trade-off 
between  localization  in  the  time  and  frequency  domains.  In  fact,  the  time  and  frequency  domain  analysis  provided  by  the 
wavelet  transform  provides  the  opportunity  to  explore  the  nature  of  transient  signals  by  representing  the  time  varying 
spectral  response  through  time-frequency  maps,  as  well  as  to  analyze  signals  for  conditions  where  responses  change 
significantly  in  amplitude  during  experiments.  Wavelets  have  found  applications  in  situations  that  require  analysis  over 
very  short  time  durations  or  where  information  is  localized,  and  have  been  combined  with  polarimetry  for 
characterization  of  biological  samples  under  transmissive  geometry.  Similarly,  fractal  analysis  is  based  on  the  original 
fractal  framework  introduced  by  Mandelbrot 15 .  Fractal  geometry  is  the  geometry  of  self-similarity  in  which  an  object 
appears  to  look  similar  at  different  scales.  The  key  concept  of  fractal  analysis  relies  on  the  fact  that  a  fractal  dimension 
can  be  considered  as  a  quantitative  measure  of  object  surface  heterogeneity  because  of  its  inherent  self-similarity 
features.  The  fractal  dimension  can  be  interpreted  as  measure  of  heterogeneity  of  a  set  of  points  on  a  plane,  or  in  space, 
for  instance,  a  measure  of  surface  roughness.  The  signals  obtained  in  this  study  are  1-D  polarimetric  signals  and  hence  a 
fractal  dimension  applicable  to  1-D  signals  must  be  applied.  There  are  several  methods  for  computation  of  fractal 
dimension  17,  such  as  the  “Box  Counting  Method,”  the  “Multiresolution  Box-counting  Method,”  the  “Katz  Method,”  the 
“Sevcik  Method,”  the  “Higuhi  Method,”  the  “Regularization  Method,”  the  “Maximum  Entropy  Method,”  and  others. 
Typically,  the  box  counting  method  is  based  on  the  method  of  covering  the  part,  whose  fractal  dimension  is  to  be 
calculated,  with  a  number  of  square  boxes  and  hence  the  number  of  boxes  N(r)  of  size  r  is  calculated  which  is  required  to 
cover  the  part.  As  the  size  of  the  square  r  (side  of  a  square)  approaches  zero,  the  area  covered  by  the  square  boxes  will 
converge  to  a  measure  of  the  curve  which  can  be  expressed  mathematically  as 

Db  =  lim(log(v(r))/log(l/r))  (1) 

In  Practice,  the  fractal  dimension  is  calculated  by  fitting  the  log-log  plot  of  log (N(r))  versus  log(l/r)  (least  squares 
method)  and  then  calculating  the  slope  of  the  plot.  The  slope  of  the  plot  gives  the  fractal  dimension  of  the  signal,  which 
gives 

log  (V  (r))  =  Z)B  log  (l/r)  +  C  (2) 

Here,  DB  is  the  fractal  dimension  or  box  counting  fractal  dimension  of  the  signal.  The  method  requires  2-D  processing  of 
the  curve  at  various  grid  sizes,  which  increases  the  computational  time.  In  this  study  we  used  an  improved  version  of  the 
box  dimension  method,  namely  the  “Multiresolution  Box-counting  Method.”17  The  applied  technique  estimates  the 
fractal  dimension  with  less  computational  time  and  is  insensitive  to  wave  amplitude,  therefore  providing  enhanced 
robustness  for  real-time  applications.  The  physical  algorithm  and  applied  metrics  of  this  study  on  the  material 
characterization  of  space  materials  is  shown  in  Fig.  1,  with  reference  to  polysilicon,  a  material  found  to  exhibit  higher 
diffused  scattering  characteristics,  depolarization,  and  fractal  dimension,  with  respect  to  amorphous  silicon. 
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Figure  1.  Physical  algorithm  and  applied  metrics 
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2.  SPACE  MATERIALS 


In  this  study,  two  space  materials  used  in  solar  panel  structures  are  examined,  namely,  amorphous  silicon  and 
poly  crystalline  silicon  (polysilicon).  Amorphous  silicon  is  an  entirely  non-crystalline  form  of  silicon  that  can  be  thought 
of  as  grains  the  size  of  individual  atoms.  On  the  other  hand,  standard  poly  silicon  consists  of  0.5  to  1  micrometer  discrete 
grains,  or  crystals,  of  mono-crystalline  silicon  that  generate  regions  of  highly  uniform  crystal  structures,  therefore,  of 
different  refractive  indices,  separated  by  grain  boundaries.  Both  amorphous  silicon  and  polysilicon  solar  cell  samples 
were  analyzed  using  confocal  microscopy.  Example  images  of  these  samples  are  shown  in  Figures  2  and  3,  respectively. 
The  amorphous  silicon  solar  sample  has  a  characteristic  size  of  approximately  2  mm,  while  the  polysilicon  solar  cell 
sample  has  grain  sizes  in  the  micrometer  range.  This  explains  the  enhanced  specular  characteristics  associated  with 
amorphous  silicon  rather  with  polysilicon,  which  instead  exhibits  pronounced  diffuse  scattering. 


Figure  2.  Confocal  microscopy  image  of  an  amorphous  silicon  solar  cell  sample 


Figure  3.  Confocal  microscopy  image  of  a  poly  silicon  solar  cell  sample 


3.  MUELLER  MATRIX  DECOMPOSITION  AS  FUNCTION  OF  THE  ASPECT  ANGLE 

Mueller  matrix  decomposition  into  a  sequence  of  three  matrix  factors,  namely,  depolarization,  retardance,  and 
diattenuation,  was  originally  proposed  by  Shih-Yau  Lu,  et  al.  32.  This  concept  can  be  expanded  for  Bidirectional 
Reflectance  Distribution  Function  (BRDF)  applications  by  expressing  them  in  terms  of  the  aspect  angle  of  an  object,  0, 
so  that 


M{e)=Mdepol{e)Mret{d)Mdiat(d) 


(3) 
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M depoi  ($)  accounts  for  the  depolarizing  effects  of  the  medium,  Mret  accounts  for  the  retardance  effects  of  linear 

birefringence  and  optical  activity,  and  Mdiat{0^  describes  the  effects  of  linear  and  circular  dichroism.  From  these 
matrices  the  diattenuation,  retardance,  and  depolarization  characteristics  of  the  medium  are  readily  determined. 

The  depolarization  is  quantified  in  terms  of  the  depolarization  index,  PD,  according  to 


PD(e)  =  Dep(M(d))  =  1 


yl(ZiJmv(0))-mZo(0) 

V3»4,  (o) 


(4) 


where  m^Q)  are  the  MM  elements  as  a  function  of  the  aspect  angle.  From  the  decomposed  retardance  matrix,  Mret  , 
the  total  retardance,  R ,  which  includes  the  effects  of  both  linear  and  circular  birefringence,  can  be  expressed  as 


/?(#)  =  COS  1 


V(M, „(<?)) 

2 

V  J 


(5) 


where  tr{^Mret  ($))  is  the  trace  of  the  retardance  matrix.  The  diattenuation,  d,  is  dependent  on  the  first  row  vector  of 

the  MM.  This  vector  describes  differential  attenuation  for  both  linear  and  circular  polarization  states  and  the 
diattenuation  can  be  expressed  in  terms  of  it  as 


d(0)  =  — f—  X  ^  m02!  (0)  +  m022  (0)  +  ml 3  Jo) 

maAO) 


(6) 


4.  WAVELET  ANALYSIS  OF  THE  POLARIZATION  STATES 

The  16  polarimetric  states  for  amorphous  silicon,  polysilicon,  Mylar,  and  Kapton  were  analyzed  using  continuous 
wavelet  transforms  (CWTs)  and  their  respective  scalograms  were  plotted.  In  this  study,  the  particular  wavelet  was 
chosen  from  the  Daubechies  family,  namely,  the  ‘db4’ wavelet.  In  order  to  understand  the  frequency  and  time  resolution 
of  the  analysis  better,  we  pick  a  line  of  coefficients  for  particular  scale  values  from  low  and  high  scale  regions 
respectively.  The  scale  values  themselves  were  picked  based  on  the  scalograms  of  the  polarization  states.  High 
frequencies  were  observed  predominantly  about  a  scale  value  of  150  while  low  frequencies  were  observed  about  a  scale 
value  of  450  on  a  range  of  5 12  scales.  In  order  to  average  noise  and  to  get  a  better  picture  of  the  signal  variation  with  the 
scale  value  chosen,  we  choose  a  range  of  scale  values  about  the  chosen  scale  line  of  interest  and  average  the  scale  line 
coefficients  in  this  range.  For  this  study,  we  choose  a  range  of  5  lines  about  the  scale  value  of  interest.  The  lowest 
frequency  was  calculated  as  follows: 


•  Each  signal  has  2500  sample  points  corresponding  to  5  cycles 

•  Each  cycle  of  signal  therefore  has  500  points 

•  Each  signal  was  captured  at  a  rate  of  105  samples/second 

•  Signal  frequency  (signal  rate/number  of  samples  per  cycle)  =  (105/500)  =  200  Hz,  which  is  also  the  modulation 
frequency  of  the  laser. 

5.  FRACTAL  DIMENSION  OF  SPACE  MATERIALS  USING  MULTIRESOLUTION  BOX 

COUNTING  METHOD 

Let  us  consider  a  continuous  time  signal.  The  signal  needs  to  be  sampled  at  a  sampling  rate  of  l/fs  at  least  twice  the 
Nyquist  rate  to  discretize  the  signal.  Let  s  =  {s(l),  s(2), . ,s(N)}  be  the  discrete  time  signal  sampled  at  a  sampling 
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frequency  of  fs  and  having  N  samples.  Each  sample  points  can  be  represented  by  (x(i),y(i)),  where  i  =  1  to  N.  The 
abscissa(x(i))  and  the  ordinate(y(i))  values  of  the  samples  are  now  known.  The  applied  algorithm  consists  of  the 
following  steps: 

Step  1:  Consider  two  discrete  samples  of  the  signal  s(i)  and  s(i+l).  The  time  difference  between  the  two  samples  is  given 
by  dt  =  x(i+l)  -  x(i),  and  the  height  difference  between  the  two  samples  is  h  =  y(i+l)  -  y(i).  The  width  of  the  rectangle 
box  required  to  cover  the  two  points  is  r  =  dt  and  the  number  of  square  boxes  of  size  r  required  to  cover  the  two  points  is 
b=|h|/r.  The  same  procedure  is  repeated  for  all  other  points  of  the  signal  with  the  value  of  h  changing  and  r  being  the 
same.  Then  the  number  of  boxes  of  size  r  required  to  cover  all  the  samples  in  the  signal  is  given  by  b(i)  =  |h(i)|/r.  The 
total  number  of  boxes  required  to  cover  the  curve  at  the  resolution  r  is  B(r)  =sum(b(i)),  where  i  =  1  to  N-l.  This  is  shown 
in  Figure  3(a). 

Step  2:  Now  the  time  resolution  is  increased  to  r  =  dt  =  x(i+2)  -  x(i),  i.e.,  the  size  of  the  box  is  increased.  The  same 
procedure  is  repeated  and  the  number  of  boxes  of  size  r  required  to  cover  all  the  samples  is  calculated. 

Step  3:  The  same  procedure  is  now  repeated  for  different  time  resolutions  and  the  number  of  boxes  required  to  cover  the 
sample  points  with  the  respective  box  size  is  calculated.  A  log-log  plot  of  B(r)  with  1/r  is  plotted.  The  plot  is  linear  fitted 
using  least  square  method  and  the  slope  of  the  linear  fit  denotes  the  fractal  dimension  of  the  1-D  signal.  The  algorithm  is 
depicted  graphically  in  Fig,  4(a)  and  4(b). 


Figure  4.  Explanation  of  multiresolution  technique  for  1-D  signals  for  two  cases:  (a)  multiresolution  technique  at  the  finest 
resolution,  and  (b)  multiresolution  technique  with  the  resolution  increased. 


To  implement  the  multiresolution  box  counting  algorithm  for  1-D  signals  from  space  materials,  the  discrete  samples  of 
the  signals  are  known  sampled  at  t  =  10"6  sec  with  N  =  2500.  Since  the  signal  is  sampled  at  very  small  time  durations,  the 
sampling  time  is  increased  to  T  =  lOt  and  thus  gives  N  =  250  samples  to  process.  The  similar  algorithm  is  applied  for  10 
grid  sizes  for  all  the  materials  (to  be  consistent)  and  the  corresponding  fractal  dimensions  of  the  signals  are  calculated. 


6.  EXPERIMENTAL  ARRANGEMENT  AND  PROCEDURES 

The  goal  of  this  study  was  to  characterize  remotely  several  types  of  space  materials  based  on  the  detection  of  distinct 
backscattered  signatures,  by  interrogating  the  objects  with  different  polarization  states  so  as  to  obtain  complete 
characterization  of  the  MM  of  the  objects  in  terms  of  depolarization,  diattenuation,  and  retardance.  The  following  space 
materials  were  studied:  a)  amorphous  silicon  (a-silicon)-based  solar  panel,  and  b)  polycrystalline  silicon  (polysilicon) 
solar  panel. 

The  experimental  arrangement  used  in  this  study  is  described  by  Figure  5. 
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Figure  5.  The  U.S.  Air  Force  Research  Laboratory  (AFRL)  multifunctional  imaging  platform 


The  U.S.  Air  Force  Research  Laboratory  (AFRL)  liquid  crystal  (LC)  polarimetric  multifunctional  imaging  platform 
comprises  a  complete  polarimeter  capable  of  deriving  the  full  16  element  Mueller  Matrix  (MM)  of  an  object  using  LC 
devices.  The  system  is  highly  automated  by  utilizing  NI  Lab  view  software  to  control  the  devices  and  measurement 
states,  and  perform  a  detailed  calibration  of  the  optical  system.  It  is  fully  reconfigurable  and  scalable  providing  enhanced 
multifunctional  surveillance  capabilities  through  dedicated  spatial,  timing  and  spectral  modules,  supported  by  advanced 
calibration  and  image  pre-processing  and  post-processing  techniques  developed  in-  house.  The  LC  polarimeter  consists 
of  the  following  components: 

1 .  Generator  arm:  The  incident  laser  light  is  initially  linearly  polarized  by  a  fixed  polarizer  at  +45°. 

2.  The  light  then  passes  through  a  compensated  variable  retarder  (polarization  rotator)  that  is  capable  of 
producing  all  linear  polarization  states,  including  but  not  limited  to  horizontal,  vertical,  and  ±45°. 

3.  The  light  then  passes  through  a  variable  retarder,  which  is  capable  of  producing  right  and  left  circularly 
polarized  light. 

4.  The  light  then  impinges  on  the  target,  and  the  diffusely  reflected  light  is  collected  by  the  analyzer  arm. 

5.  Analyzer  arm:  The  analyzer  arm  consists  of  the  same  components  as  the  input  arm  except  in  reverse  order  and 
with  the  fixed  polarizer  set  at  -45°.  The  experimental  arrangement  of  the  generator-analyzer  arms  is  shown  in 
Figure  6. 


Figure  6.  Experimental  geometry  of  the  backscattering  polarimetric  system 
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In  these  sets  of  experiments,  1 065  nm  laser  pulses  with  a  PRR  of  200  Hz  were  transmitted  through  the  generator  system 
illuminating  the  space  object  material.  The  light  was  then  backscattered  in  the  direction  of  the  analyzer  arm  and  detected 
by  a  New  Focus  2151  femtowatt  photodetector.  The  waveforms  were  acquired  through  a  NI  Labview  data  acquisition 
system  (DAS),  and  then  processed  using  Excel  and  Matlab  subroutines.  Amorphous  silicon  and  polysilicon  solar  panels 
were  mounted  on  a  rotating  stage.  As  shown  in  Figure  6,  the  receiver  was  placed  off-axis  with  respect  to  the  normal 
incidence  of  laser  beam  on  the  solar  panel  at  a  fixed  scattering  angle  of  2.46°.  Rotation  was  achieved  by  way  of  a  stepper 
motor  with  a  resolution  of  0.01°. 

6.1  Calibration  Procedure 

A  technique  called  the  Method  of  Swings  [35]  was  used  in  conjunction  with  the  Null  Intensity  Method  [19].  This 
involves  setting  Polarizer( Analyzer)  at  some  angle  +/-  d0  about  the  suspected  null.  When  the  intensities  are  identical  the 
true  location  of  the  Null  is  known.  The  full  calibration  procedure  applied  to  this  study  is  the  following: 

1)  Use  Null  Intensity  Method  to  determine  the  LC  voltage  which  corresponds  to  the  desired  rotation  or  retardation 
to  within  1/1 0th  of  a  volt. 

2)  Use  the  “Method  of  Swings”  to  further  adjust  the  LC  voltage  for  maximum  precision. 

3)  Verify  the  proper  calibration  of  the  state  by  varying  the  0  of  Polarizer(A)  according  to  the  method  described  in 
[19]. 

A  comparison  of  the  experimentally  estimated  1 6  polarization  states  with  reported  accuracy  (for  the  calculation  of  MM 
of  Air,  Linear  Horizontal  Polarizer  (LHP)  and  Linear  Vertical  Polarizer  (LVP)  with  ideal  and  measured  values)  is  shown 
in  Table  I. 


Table  I.  Calibration  and  System  Accuracy 


Ideal  Mueller  Matrix 

Measured  Mueller  Matrix 

Element  Vari  a  ntes 

Element 

Air 

LVP 

LHP 

Air 

LVP 

LHP 

Air 

LVP 

LHP 

mil 

1 

1 

1 

1.000 

1.000 

3  1.000  T 

0.0 

0.0 

O.O 

m!2 

0 

*1 

1 

0.001 

-0.969 

0.974 

1.3E-6 

9S5.4E-6 

666.7E-6 

ml} 

0 

0 

0 

-0.043 

0.057 

-0.029 

l.SE-3 

3.3E-3 

S3S.7E-6 

mI4 

0 

0 

0 

-0.012 

-0.001 

-0.011 

141. 3  E- 6 

760.3E-9 

110.8  E-6 

m21 

0 

'1 

1 

0.001 

*0,968 

0,971 

1-1E-6 

l.QE-3 

813.5E-6 

111  22 

1 

1 

1 

0.992 

0.960 

0.971 

60.0  E-6 

1.6  E-3 

827.2E-6 

m2J 

0 

0 

0 

0.002 

-0.058 

■0,028 

2.8  E-6 

3.3E-3 

771.1E-6 

m24 

0 

0 

0 

0.045 

0.001 

-O.OW 

2.0E-3 

258.7E-9 

93.2E-6 

mil 

0 

0 

0 

0.070 

-0.096 

0.011 

4.9  E- 3 

9.2E-3 

115 .9  E-6 

m32 

0 

0 

0 

-0.016 

0.094 

0,011 

269. 8  E-6 

8.8  E-3 

117. 6  E-6 

mil 

1 

0 

0 

0.991 

0  019 

0.005 

S6.4E-6 

343.3E-6 

24.7E-6 

mJ4 

O 

0 

0' 

-0.053 

0.015 

-0.001 

2-8  E- 3 

239.5E-6 

1.2E-6 

m4l 

0 

0 

0 

0.039 

0.001 

-0.009 

1.5  E-3 

1.0  E-6 

82.8  E-6 

m42 

O 

0 

0 

0.009 

-0.002 

-0009 

82.5  E-6 

6.1E-6 

81.3E-6 

m43 

0 

0 

0 

-0.022 

0.023 

-0.002 

500.2  E-6 

533.9E-6 

4. 8  E-6 

m44 

1 

0 

0 

0.989 

0.023 

0.003 

123. 5  E-6 

527.2E-6 

10.4E-6 

STANDARD  DEVIATION: 

29.9  E-3 

43.2E-3 

16. 9  E-3 

7.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Experiments  were  performed  under  quasi-monostatic  backscattered  polarimetric  geometry  using  amorphous  silicon  and 
polysilicon  material  samples.  Backscattered  polarimetric  intensities  as  a  function  of  the  aspect  angles  of  the  solar  panel 
samples  were  recorded  and  the  Mueller  matrices  of  the  samples  were  estimated.  By  applying  the  MM  decomposition 
technique  (Eq.  3),  the  depolarization,  diattenuation,  and  retardance,  matrices  for  a-silicon  and  polysilicon  were  estimated 
and  plotted  as  function  of  the  aspect  angle,  as  shown  in  Figures  7(a)  and  7(b),  respectively.  The  experimental  results  of 
Fig.  7  indicate  that  polysilicon  exhibits  higher  depolarization  than  amorphous  silicon.  This  is  attributed  to  its  pronounced 
diffuse  scattering  characteristics  because  of  its  microcrystalline  structure,  presence  of  defects  and  none  uniformity.  On 
the  other  hand,  the  retardance  exhibits  a  Lambertian  pattern  in  both  samples.  The  diattenuation  of  the  amorphous  silicon 
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and  polysilicon  increases  initially  with  increasing  the  aspect  angle  of  the  object,  while  polysilicon  exhibits  a  distinct 
pattern  which  was  confirmed  through  repetitive  measurements;  the  symmetric  decrease  of  the  diattenuation  at  larger 
aspect  angle  of  the  polysilicon  may  be  attributed  to  the  drop  of  the  signal-to-noise  ratio  because  of  polysilicon’s  intrinsic 
pronounced  scattering  characteristics. 


Figure  7.  Decomposition  matrices  of  two  types  of  space  materials:  (a)  amorphous  silicon,  and  (b)  polysilicon 


Analysis  of  the  experimentally  determined  polarization  states  of  the  space  materials  was  performed  using  wavelet 
transforms  according  to  [12].  Scalograms  were  calculated  for  an  aspect  angle  ranging  between  -0.07  to  +0.07  degrees, 
with  scales  from  2  to  512,  in  step  size  of  2;  representative  examples  for  selected  cases  are  shown  in  Figure  8. 


(c)  (d) 

Figure  8.  Scalograms  corresponding  to  amorphous  silicon  and  polysilicon  at  an  aspect  angle  of  -0.07  degrees.  Amorphous 
silicon:  (a)  copolarized  geometry  and  (b)  crosspolarized  geometry.  Polysilicon:  (c)  copolarized  geometry  and  (d) 
crosspolarized  geometry 


Distinct  signatures  are  observed  reflecting  mainly  the  specular  nature  of  amorphous  silicon  and  the  diffuse  scattering 
characteristics  of  polysilicon;  diffuse  scattering  contributes  to  blur  or  “ghost”  patterns,  at  the  low  and  high  frequencies.  It 
has  been  observed  that  blur  or  “ghost  patterns  increases  with  increasing  the  aspect  angle  of  the  object,  as  diffuse 
reflectance  becomes  dominant;  again,  polysilicon  exhibits  these  effects  to  a  greater  degree  due  to  its  higher 
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depolarization  and  scattering  yield.  In  search  of  distinct  signatures  from  these  two  space  panel  materials,  a  line  of 
coefficients  for  particular  scale  values  from  low-scale  (high-frequency)  and  high-scale  (low-frequency)  regions, 
respectively,  was  selected  and  plotted.  The  scale  values  themselves  were  picked  based  on  the  scalograms  of  the 
polarization  states.  High  frequencies  were  observed  predominantly  about  a  scale  value  of  150  while  low  frequencies 
were  observed  about  a  scale  value  of  450  on  a  range  of  512  scales.  In  order  to  average  noise  and  obtain  a  better 
representation  of  the  signal  variation  with  the  scale  value  chosen,  a  range  of  scale  values  about  the  chosen  scale  line  of 
interest  (148  to  152)  or  (448  to  452)  was  selected;  specifically,  a  range  of  5  lines  about  the  scale  value  of  interest  was 
chosen.  An  average  of  high-scale  (low-frequency)  line  coefficients  for  amorphous  silicon  and  polysilicon,  both  for 
selected  polarization  stated,  is  shown  in  Figure  9. 


(a)  (b) 

Figure  9.  Average  of  high-scale  (low-frequency)  line  coefficients  for  amorphous  silicon  and  polysilicon:  (a)  amorphous 
silicon;  HH,  HV,  PP,  RH,  RR  states,  and  (b)  polysilicon:  HH,  HV,  PP,  RH,  RR  states 


Signal  similarity  has  been  quantified  in  terms  of  the  correlation  coefficients  of  different  polarimetric  scalograms  each 
one  representing  one  of  the  16  polarization  states  of  a  space  material  sample  at  a  fixed  aspect  angle.  A  pictorial 
representation  of  the  autocorrelation  between  polarization  for  amorphous  silicon  and  polysilicon-based  solar  panel 
samples,  at  two  different  aspect  angles,  -0.07  and  0  degrees  are  shown  in  Fig.  10. 


(a)  (b)  (c)  (d) 

Figure  10.  Autocorrelations  among  signals  representing  the  16  states  for  amorphous  silicon  and  polysilicon  at  two  different 
aspect  angles:  (a)  amorphous  silicon  at  0  degrees,  (b)  polysilicon  at  0  degrees,  (c)  amorphous  silicon  at  -0.07  degrees,  and 
(d)  Poly  silicon  at  -0.07  degrees 
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The  x  and  y  axes  represent  the  16  polarization  states  linearly  and  uniformly  distributed.  The  white  color  represents  a 
value  of  1  indicating  a  100%  likeness,  while  progressively  darker  shades  indicate  a  proportionally  reduced  degree  of 
likeness.  In  other  words,  the  whiter  regions  represent  high  correlations.  Each  pixel  value  represents  the  correlation 
coefficient  between  polarization  states  at  fixed  material,  and  fixed  angle.  By  observing  Figure  9,  it  is  seen  that 
amorphous  silicon  and  poly  silicon  exhibit  more  signal  similarity  (whiter  regions)  at  0  degrees  than  at  -0.07,  within 
certain  regions  of  interest  (ROI).  This  would  be  explained  based  on  diffused  scattering  characteristics  of  poly  silicon 
which  becomes  pronounced  with  increasing  aspect  angles.  The  fractal  dimension  is  calculated  by  fitting  (best  linear  fit) 
the  log-log  plot  of  log(N(r))  and  log(l/r)  (least  squares  method)  and  then  calculating  the  slope  of  the  plot.  The  slope  of 
the  plot  gives  the  fractal  dimension  of  the  signal,  according  to  Eq.  2.  The  logarithmic  plots  of  box  counts  (N(r))  versus 
resolution  1/r  for  amorphous  silicon  and  polysilicon,  are  shown  in  Figure  11. 

(a).  Amorphous  Silicon  at  0.14  degrees:  Backscattered  signal  and  Box  counts  vs.  resolution 

Backscattered  signal(HH)  of  Amorphous  Silicon  angled  at  0.14  degrees  Logarithmic  Plot  of  BoxCounts  Vs  resolution 


(b)  Polysilicon  at  0.14  degrees:  Backscattered  signal  and  Box  counts  vs.  resolution 


Backscattered  signal(HH)  of  Polysilicon  angled  at  0.14  degrees 


Logarithmic  Plot  of  BoxCounts  Vs  resolution 


Figure  1 1.  (a)  Backscattered  contributions  and  logarithmic  plots  of  box  counts  (N(r))  versus  resolution  1/r  at  an  aspect  angle 
of  0.14°  for  amorphous  silicon  under  collinearly  polarized  geometry;  (b)  backscattered  contributions  and  logarithmic  plots 
of  box  counts  (N(r))  versus  resolution  1/r  at  an  aspect  angle  of  0.14°  for  polysilicon  under  collinearly  polarized  geometry 

A  comparison  of  fractal  dimensions  among  different  space  material  for  selected  polarization  states  and  aspect  angles  is  shown  in 
Table  II: 
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Tablell  Comparison  of  fractal  dimensions  among  different  space  materials 


Materials/ Angles 

0  degrees 

0.14  degrees  (predominantly  diffused 
reflectance) 

HH 

HV 

HH 

HV 

Amorphous  Silicon 

1.189644534 

1.435106471 

1.289284645 

1.427184698 

Mylar 

1.203158205 

1.222856222 

1.197328496 

1.269531896 

Kapton 

1.224981172 

1.355107509 

1.290935866 

1 .389688974 

Polysilicon 

U48899471 

1.453745597 

1.335139104 

1.512163558 

The  fractal  dimensions  of  the  amorphous  silicon,  polysilicon,  Kapton,  and  Mylar,  are  plotted  for  distinct  polarization 
states  in  Figure  12,  at  two  different  aspect  angles. 


States  at  two  different  angles 


■  Amorphous  Silicon 

■  Mylar 

■  Kapton 

■  Polysilicon 


Figure  12.  Plot  of  fractal  dimension  of  space  materials  at  two  different  polarimetric  states  and  two  different  aspect  angles 


8.  CONCLUSION 

The  experimental  results  of  this  study  indicate  that  media  with  highly  diffuse  scattering  characteristics,  such  as 
polysilicon,  exhibit  higher  depolarization,  and  fractal  dimension,  with  respect  to  amorphous  silicon  which  exhibit  rather 
pronounced  specular  characteristics.  The  outcome  of  this  experimental  study  indicates  that  polarimetric  wavelet  fractal 
principles  applied  to  ladar  return  signals  may  provide  an  enhanced  capability  to  characterize  and  discriminate  unresolved 
space  objects. 
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